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Abstract 

The effect of flame stretch, composed of strain and curvature, plays a major role in the propagation of 
turbulent premixed flames. Although all forms of stretch (positive and negative) are present in turbulent 
conditions, little research has been focused on the stretch due to curvature. The present study 
the Markstein number (which characterizes the sensitivity of the flame propagation speed to the impo 
stretch rate) for an inwardly-propagating flame (IFF). This flame is of interest b-" is negatively 
stretched, and is subjected to curvature effects alone, without the competing effec s _ 

extension of our previous work, microgravity experiments were run using a vortex-flame ‘^racUon to 
create a pocket of reactants surrounded by an IPF. Computations using the RUN-1DL code of Rogg 
were also performed in order to explain the measurements. . . 

It was found that the Markstein number of an inwardly-propagating flame, for both the microgravity 
experiment and the computations, is significantly larger than that of an outwardly-propagatmg flame. 
Further insight was gained by running the computations for the simplified (hypothetical) cases of one- 
step chemistry, unity Lewis number, and negligible heat release. Results provide additional evidence that 
the Markstein numbers associated with strain and curvature have different values. 
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INTRODUCTION 


Because turbulence exerts both curvature and strain on a premised d.,- .kjlu hU 

studies of how flames react to both of these conditions [1-11]. However, most of these u. L i«-s 1 m\c Ibru-od 
on strained flames, such as the counterfiow flame case [ij or flames with positive strain and curvature such 
as the spherical outwardly-propagating flame (OPF) [2-5]. The present work instead focuses on one type 
of negatively-stretched flame, namely the inwardly- propagating flame (IPF). The IPF has the desirable 
property that the stretch rate (normally composed of both strain and curvature) consists of curvature effects 
alone. The complicating effects of strain are not present, since the gas ahead of the IPF is at rest. Ideally, 
it would be useful to quantify the Markstein number, which represents the sensitivity of the propagation 
speed to the imposed stretch rate, for four conditions: positive strain alone (the counterflow case), positive 
curvature alone, negative strain alone, and negative curvature alone (which is the IPF condition). Simple 
theory predicts that the Markstein number, A/a, will be the same for all four of these conditions (at the 
same fuel-air equivalence ratio 4>). However, results of Bradley et al. [5], Echekki and Mungal [6], Poinsot 
et al. [7], Durox et al. [8], and Ibarreta and Driscoll [9] indicate that there may be different values of Ma 
associated with each of these four cases. 

The present work builds upon a previous experimental study [9] conducted at 1-g, in which it was shown 
that Ma for an IPF is two to three times larger than that of a corresponding OPF. In the present work the 
experiment is operated in microgravity, which removes the complications of buoyancy and flame-produced 
vorticity [12, 13], and a full set of numerical computations are conducted, using the RUN-1DL code of 
Rogg [14, 15] with complex chemistry, to better explain the experimental results. Computed values of Ma 
are reported for both the OPF and the IPF cases, and are compared to the measurements. An advantage 
of the experiment is that a nearly spherical IPF is formed in the latter stages of the burnup of a pocket 
of reactants. This spherical geometry allows for accurate measurement of flame speed since the reactants 
are essentially at rest. Previous measurements of negatively-stretched flames have been reported for more 
complicated geometries [6-8, 10, 11, 16]. 

The Markstein number (Ma) that was measured relates the stretched burning velocity, S u , to its un- 
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stretched value, 5 t f , through the following relation: 


This equal ion is based or* l)v' .•.rkduu] •dys;^ of \ T ■ -k'-f • ' ; T i - J T); h - *' ■ Vfi i lur « :-f K-'flovi** 

number: 

Ka — K/ [( S u ) 2 /£>„] (2) 

Where /C is the stretch rate; Kwon et al. [2] have suggested the normalization factor that appears in Eq. 2, 
which depends on the local burning velocity of the stretched flame ( S u ). D u is the binary diffusivity of 
the fuel into nitrogen. The burning velocity is defined as the speed of flame propagation relative to the 
unburned reactants. The stretch rate [17] equals A~ l dA/dt\, where A is an infinitesimal area of flame 
surface. Note that other researchers have chosen to nondimensionalize the stretch rate using the unstretched 
velocity, 5,?, rather than S u . This leads to a relation that is slightly different from Eq. 1, for which the 
right side is (1 — Ma Ka). Both relations yield exactly the same values of Ma in the limit as Ka approaches 
zero. However, experiments show that Eq. 1 provides a good fit to experimental data (as shown below) for 
moderate values of Ka y and thus provides an accurate method to determine values of Ma. 

MICROGRAVITY EXPERIMENT 

The experimental setup is similar to a previous 1-g experiment [9] that was modified to be operated at 
microgravity conditions in the NASA Glenn 2.2 Second Drop Tower. A premixed propane-air flame, at an 
equivalence ratio (<£) of 0.54, is interacted with a vortex propagating in the opposite direction, causing the 
flame to wrap around the vortex. The distorted flame surface eventually pinches off from the main surface, 
creating a pocket of reactants surrounded by an IPF. Long after the vortex has been convected away and 
has decayed, the IPF burns from an initial diameter of approximately 10 mm to a final observed diameter 
of less than 1 mm. The combustion chamber is 11.4 cm by 11.4 cm by 23.5 cm and has a loudspeaker at one 
end to force gas through a 2.54 cm orifice, creating the toroidal vortex. In a series of events programmed 
into a Tattletale controller, an electrical connection is broken which signals the beginning of the drop; after 
a 500 ms wait time, exhaust ports are opened to prevent pressure buildup and the flame is ignited using a 
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spark. After 200 ms the speaker is pulsed and the pocket formation occurs 300 - 400 ms later. Final pocket 



a Kodak RO digital camera. Figure 1 shows the digitized image of the flame boundary at four different times. 
Since the actual flame is not perfectly spherical, the boundary at each frame was fitted with a polynomial 
and then divided into 400 flame segments to determine the local normal vector and the burning velocity 
( S u ), which is: S u — AX n / At. The distance from one segment of the flame front to the flamefront in the 
next image is AX n , as measured along the normal direction, and A t is the time between images. Candel 
and Poinsot [18] have shown that the stretch rate ( K ) is -n • (n • V)u + V • u + S u / R c . For the IPF the 
velocity of the gas (u) inside the pocket is zero due to symmetry, so the stretch rate of an IPF is: 

K = S u /R c , (3) 

The radius of curvature, R c , was calculated for each segment from its in-plane and out-of-plane compo- 
nents, assuming that the pocket is symmetrical about the vertical axis. The measured Markstein number 
was determined for the IPF by combining Eqs. 1-3 and solving for S u , which yields: 

S u / = 1 - Ma 5 d / R c (4) 

The unstretched flame thickness is defined as D u /S£. This gives us a simple relation between the 
burning velocity and the flame curvature (normalized by <5 D °). Note that R c is a negative number for the 
IPF geometry considered. 

NUMERICAL SIMULATIONS 

To explain the experimental findings, numerical simulations of both inwardly and outwardly propagating 
spherical flames (with complex chemistry) were generated using the RUN-1DL code of Rogg [14, 15]. This 
computer code solves the time-dependant conservation equations in spherical coordinates. An adaptive grid 
tracks the flame front as it moves in time. Experimental conditions were matched by considering a lean 
propane-air flame with initial conditions of 300 K and 1 atm. Symmetric boundary conditions are enforced 
both at the origin and at the edge of the domain, so that the gradients of temperature and of species mass 
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fractions at the boundaries are set to equal zero. The propane-air detailed chemical mechanism of Djavdan et 

ai. [19] was used which includes 29 specie*? and 1 rmr! ^ ; 

agrees with experimental data of Egoifopoulos et ai. [20] to within 20%. To examine if the same trends in the 

using the detailed chemistry of Smooke et al. [21] that includes 16 species and 46 reactions and predicts an 
unstretched burning velocity that agrees with experimental data [22] to within 10%. In addition, the role 
of intermediate species was investigated by removing intermediates from the calculation entirely; values of 
Ma were calculated for simple-one step methane-air chemistry and were compared to the complex chemistry 
results. For the one-step chemistry, an activation energy of 50 keal/mol and an Arrhenius pre-exponential 
factor of 4.4 ■ 10“ 15 yield a stoichiometric unstretched burning velocity that agrees with experiments [22] to 
within 10%. 

The computed flame boundary was defined as the radial location of maximum temperature gradient, 
since this is comparable to the boundary defined by the shadowgraph images obtained in the experiment. 
For the IPF computations, the stretch rate was determined using Eq. 3 and the local burning velocity was 
determined using: 

Su = d V u jdt (5) 

since the gas velocity ahead of the IPF is zero; T u is defined as the radial position of the flame surface. The 
computed Markstein numbers for the IPF were determined by combining Eqs. 4 and 5. 

For the computed outwardly-propagating flame, the gas velocity ahead of the flame is not zero, so Eqs. 
3-5 are not applicable. Instead, the stretch rate and the local burning velocity for the spherical OPF can 
be shown [2] to be [(2/ T u ) • dV u jdt} and [(p&/pu) - d T u /dt\, respectively. Once K and S u are computed 
in this manner, the Markstein number of the OPF is found using Eq. 1 and 2. 

Several researchers have published corrections for the Markstein number of the OPF related to accu- 
mulation and traverse fluxes along the flamefront due to the finite flame thickness [1,5,23,24]. There is 
disagreement, however, between each of the corrections which depend heavily on the definition of the flame 
thickness and therefore add ambiguity to the results. It was decided not to apply the corrections in the 
current work, so that the flame is assumed to be infinitely thin with respect to the radius of curvature. It is 
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important to note that none of the corrections apply to the IPF configuration. 


RESULT'S 

Figure 2 shows some results of the microgravity IPF experiment, for which the propane-air equivalence 
ratio was 0.54. The slope of the curve in Fig. 2 is (—Mi), in accordance with Eq. 4. Figure 2 was obtained 
by compiling about 2000 data points from one drop; the flame boundaries were divided into segments and 
the local curvature and burning velocity were calculated as described earlier. The data were then grouped 
according to curvature (1 / R c ) and plotted. Each solid symbol represents the mean burning velocity for the 
points in each of the curvature groups. The error bars represent the range from -hi to -1 standard deviations 
(a) from the mean of the group. Note that the uncertainty in determining the mean value (i.e. the solid 
symbol) is much less than the height of the error bars drawn. Figure 2 shows that, despite the large amount 
of scatter associated with this difficult experiment, there appears to be a linear relation between the variables 
plotted; therefore Eq. 4 is valid. The slope of the data in Fig. 2 yields a Markstein number of 12.1. 

Figure 3 shows some numerical results obtained using the RUN-1DL for both OPF and IPF cases, with 
complex propane-air chemistry at $ = 0.6. The slope of the curve for the IPF (i.e. negative Karlovitz 
numbers) yields a Markstein number of 11.5, which is within the range of values that were measured in the 
microgravity experiment. The difference in slopes seen in Figure 3 also indicate that the Markstein number 
of the computed IPF is larger than that of the corresponding OPF. 

Comparisons are made in Fig. 4 between the Markstein numbers obtained in the experiment, the com- 
puted IPF, and the computed OPF cases. There is satisfactory agreement between the experimental IPF 
values (solid diamonds) and the computed IPF values (solid triangles). Also shown are the results of Sun et 
al. [1], represented by the solid circles. Sun et al. also performed IPF computations for complex propane-air 
chemistry, and their results lie within 20% of the present values. The OPF results of the present work are 
shown as open triangles. Also shown are OPF results from experiments [3,4,25,26] and from computations. 
It can be concluded from Fig. 4 that in all cases the Markstein numbers associated with the IPF are 2-3 
times larger than those for the OPF. 

The numerical results are useful in that they show how positive and negative stretch affect the profiles 
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of species concentrations and temperature across the flame. Two examples are given by the IPF profiles 

that ported f*« com pared to (Vm nn«+ro*chod flame). This trend is expected: the value of Ma is positive 

and the negative stretch exerted causes an increase in the burning velocity, as given by Eq. 1. A similar and 
expected increase in the OH and H radical concentrations is also observed. 

In Fig. 6 the reverse trend occurs, as expected. As positive stretch is exerted, the temperatures within 
the stretched case (solid line) now lie below those of the unstretched case (dotted line). This decrease in 
flame temperature is tied to a decrease of radical concentrations such as OH and H, and also corresponds to 
a decrease in burning velocity. 

An interesting difference between the IPF and OPF is that there is a buildup of hydrogen molecules 
ahead of the IPF (see Fig. 5), but no such buildup occurs for the OPF. This effect is due to the zero-gradient 
boundary condition at the origin, which allows the H 2 diffusing upstream of the IPF flame to accumulate. 
This mimics the flame-flame interaction that would occur in the IPF during the very late stages of pocket 
burnup. The influence of this effect was quantified by changing the boundary condition and setting the mole 
fractions constant at r = 0. No significant change in the Markstein number was observed, thus flame- flame 
interaction was assumed to be trivial for the range of stretch values we are interested in. 

To determine if the above results are unique to propane-air flames, or are also observed for other hy- 
drocarbon fuels, the computations were repeated for methane-air flames with complex chemistry. Figure 7 
shows that the methane-air flames are similar to propane-air flames in that the values of Ma for the IPF 
(solid triangles) are more than 2-3 times larger than for the OPF. The values of Ma computed for the OPF 
cases are in agreement with the OPF values measured by Tseng et al. [3,25], which are represented by the 
open circles in Fig. 7. 

Since the computational results agree with the experimental finding that values of Ma are larger for the 
IPF than for the OPF case, it was decided to determine if radicals and intermediates cause this difference. 
Therefore radicals and intermediates were removed by running the RUN-1DL code using one-step chemistry. 
The one-step chemistry results are compared to the complex chemistry results in Fig. 7. For this comparison, 
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methane-air mixtures were considered, since a number of researchers have proposed global or quasi-global 

reaction schemes for rv**« rX fu?-’ r .vw ■ n*. .... : 

luis computation. Figure 7 shows that one-step chemistry does predict values of Markstein number (the solid 
vlM.mou'-.i sym: i. hat are ill agree:, ;ern. AV.h. die coniplex i;i > ..uU?.rv o..s;xIis (de- Li 

ihe one-step chemistry results also predict that Ma for the IFF case is larger than that for the OFF case. 
It is concluded that the presence of radicals and intermediate species is not the primary factor that causes 
the difference in Markstein number of the IPF and OPF cases. 

Another computation was performed to determine the role, if any, of the Lewis number of the deficient 
reactant. The Lewis number ( Le ) is defined for each species as the ratio of the thermal diffusivity to the 
mass diffusivity of the species. The RUN-1 DL code was run for a simple one-step methane-air chemistry 
using values found by Coffee [28]; the Lewis number of each species, except for the fuel, was set to one. The 
fuel Lewis number (Le F ) was varied to observe its effects on the Markstein number. Figure 8a shows that 
for these simplified conditions, the Markstein number is, as expected, directly related to the Lewis number 
of the deficient reactant. There is still, however, an unexplained difference in the IPF and OPF Markstein 
numbers. This difference increases for larger Lewis numbers. 

Finally, the role that thermal expansion plays in determining the Markstein number was investigated. 
For an OPF, thermal expansion caused by heat release increases the velocity of the reactants ahead of the 
flame and creates a strain field, since the expanding flame acts as a piston. To study the flame without this 
complicated effect, the heat release rates for both IPF and OPF cases were set to be nearly zero. RUN-1DL 
computations were obtained with single-step chemistry and for sufficiently small enthalpy of combustion 
such that the temperature rise across the flame was only 10%. Other constants, such as activation energy, 
were adjusted to yield an unstretched burning velocity similar to that of a methane-air mixture at the same 
equivalence ratio (0.66). Results for this cool flame simulation are seen in Fig. 8b. These results show that 
Markstein number is identical for both the IPF and OPF cases over a wide range of conditions if heat release 
is absent. Therefore the differences in Markstein number between the OPF and IPF cases can be attributed 
to different values of Markstein number associated with the strain field and the curvature effects. Other 
studies which have investigated these differences have been conducted by Poinsot et al. [7] and Bradley et 
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Experimental u nr ertainties for this study have been documented in Ref. 9. The computational uncer- 
tainties were determined by varying the number of grid points, the initial conditions, and the domain size. 
The number of radial grid points was chosen to be between 100 and 140; a sample calculation was performed 
with half the normal number of grid points and the resulting value of Ma changed by less than 5%. To spec- 
ify the initial conditions, the entire domain is initially composed of only reactants, and a Gaussian-shaped 
temperature profile is imposed at outermost radial location of the domain, causing the mixture to ignite and 
the flame to burn inward. The excess enthalpy associated with this ignition source rapidly diffuses radially 
out of the domain. To insure that the results are independent of the initial temperature profile chosen, the 
value of the maximum initial temperature was raised from 1700 K to 2400 K. For both of these values the 
flame profiles quickly adjusted to become identical, and the resulting values of Ma differed by less than 1%. 
The domain size ideally would be infinite, but it was found that values of Ma varied by less than 3% if the 
domain size was varied from 5 cm to 9 cm, so values in this range were used. 

CONCLUSIONS 

1. Markstein numbers of a inwardly-propagating flame (IPF) were measured in microgravity using a 
unique experiment. This flame is of interest because it is subjected to negative stretch rates, and because 
only the curvature effect is present; the complications associated with strain are absent. The measured values 
of Ma for the IPF are significantly larger than previous measurements reported for the OPF case. 

2. Computed Markstein numbers also were obtained using the RUN-1DL code and complex chemistry. 
Results agreed with the experiment. Computed values of Ma for the IPF were significantly larger than 
computed for the OPF case, in agreement with the measurements. 

3. When the thermal expansion is artificially removed from the computations, both IPF and OPF cases 
are subjected to curvature effects alone, and strain effects are absent; the resulting computed Markstein 
numbers are identical. This indicates that the Markstein numbers associated with strain and curvature 



differ, which has been suggested previously in worV bv Poinsot et ah [71- B radio v ot, ah W and others. 


1,1 T 5 * ;i ' hfFusvai of the various species, based on computations for which these quantities 

were varied. 
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Figure 1: Experimental Same contours of an inwardly propagating flame (IFF) at microgravity conditions 
are plotted tor sever ai tunes; time step is 4 ms. The boundary is first determined using a ;o > . , * 

a'i'- du'n fhr.pM with polynomials The flame q>eed ^ rui* -u ]a?eo ].»> allv esi ■ * ’ ■ 
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Figure 2: Experimental Results (Propane - $ = 0.54): Compilation scatter plot of non-dimensional 
flame speed ( S u /S ° ) versus non-dimensional curvature (S° n /R c ). The size of the circles indicate the number 
of points averaged at that curvature location. Tiie error bars extend one standard deviation above and below 
the mean for the points averaged. The slope of the least-squares fit line through point (0,1) is equal to -A/a, 
as shown in Eq. 4. 


Figure 3: Numerical Results (Propane - = 0.6): Computed [S°/S u ) versus non-dimensional stretch 

rate (A a) for OPFs (□) and IPFs (a). The Markstein number is defined as the slope of the best linear fit 
(shown as a solid line) from Eq. 1. Notice that the IPF Markstein number is much larger, such that the 
slopes clearly do not match. 


Figure 4: Propane- Air Results: Compilation plot showing numerical and experimental results obtained 
for both IPF (solid symbols) and OPF (open symbols) flames. Experimental microgravity results of current 
study are shown (♦). Current numerical results for IPFs and OPFs are shown as triangles (A and A ). 
This set of data is compared to numerical data obtained for IPF and OPF flames of Sun et al. [1], shown 
as circles (o and • ). Note that their OPF data were modified by using our definition of burning velocity. 
Other experimental results for OPFs of Tseng at al. [3], Taylor [26] and Palm-Leis et al. [4] are also shown 
for comparison (0). All available data show a definite difference between the Markstein values associated 
with OPFs and IPFs. 


Figure 5: Propane- Air IPF Case (4> = 0.6): Computed temperature and molar concentrations of choice 
species across the outwardly propagating flame. The reference unstretched case is shown as the dotted lines. 
The solid lines represent the flame structure for the stretched case (K - -704 s~ 1 ). The x = 0 location 
corresponds to the point of maximum heat release. The origin is located to the left of the plot. Notice that 
stretch causes a moderate increase in final temperature and intermediate species such as H, OH and CO; 
this increase is linked to a large increase in burning velocity (because of the positive A/a). The increase in 
intermediate species is also associated with an increase in heat release. For the stretched case there is also a 
very noticeable increased concentration of H 2 ahead of the flame. The flame thickness, as defined from the 
temperature profile, is greatly diminished for the stretched flame. 


Figure 6: Propane- Air OPF Case (4> — 0.6): Computed temperature and molar concentrations of choice 
species across the outwardly propagating flame. The reference unstretched case is shown as the dotted lines. 
The solid lines represent the flame structure for the stretched case (K = +113s -1 ). The x = 0 location 
corresponds to the point of maximum heat release. The origin is now located to the right of the plot. Notice 
that stretch causes a small decrease in final temperature and intermediate species; this decrease is linked to 
a decrease in burning velocity (because of the positive A/a). The decrease in intermediate species is also 
associated with an decrease in heat release. 
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Figure 7: Methane-Air Results: Compilation plot showing numerical results obtained for both IPF (solid 
symbols) and OPF (open symbols) flames. Numerical results for complex chemistry IPFs and OPFs are 
shown as triangles (A and A ). Results obtained for simple step chemistry are shown as diamonds ( 0 and 
♦ ). Note that IPF Markstein values are larger then OPF values. Experimental OPF data of Tseng et 
al. [3,25] are shown for comparison ( o ). 


Figure 8: Lewis number dependence for simple chemistry ( Methane-air; $ = 0.66 ): Numerical 
results of Markstein numbers versus Lewis number of the fuel for methane/air flames using simple one- 
step chemistry. Both plots show results for IPFs (♦) and OPFs (0). Figure 8(a) shows the dependence 
of Markstein number on the fuel Lewis number for the simple mechanism. Notice there is still a large 
discrepancy between values obtained for IPFs and OPFs. Figure 8(b) shows results when heat release is 
lowered to negligible values. In this case, both flame configurations are purely curved. Both Markstein values 
collapse onto each other for a wide range of Lewis numbers. 
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X-Coordinate (cm) 


Figure 1: Experimental flame contours of an inwardly propagating flame (IPF) at microgravity conditions 
are plotted for several times; time step is 4 ms. The boundary is first determined using a shadowgraph system 
and then fitted with polynomials. The flame speed is calculated locally using the normal distance between 
surfaces, AX n) and the time difference, At. Only 1 of every 2 surfaces obtained is shown for clarity. 
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Figure 2: Experimental Results (Propane - $ = 0.54): Compilation scatter plot of non-dimensional 
flame speed (S u /S°) versus non-dimensional curvature (6° D /R C ). The size of the circles indicate the number 
of points averaged at that curvature location. The error bars extend one standard deviation above and below 
the mean for the points averaged. The slope of the least-squares fit line through point (0,1) is equal to -Ma, 
as shown in Eq. 4. 
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Figme 3: Numerical Results (Propane - $ = 0.6): Computed (S°/S u ) versus non-dimensional stretch 
rate (Ka) for OPFs (□) and IPFs (a). The Markstein number is defined as the slope of the best linear fit 
(shown as a solid line) from Eq. 1. Notice that the IPF Markstein number is much larger, such that the 
slopes clearly do not match. 
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Figure 4: Propane- Air Results: Compilation plot showing numerical and experimental results obtained 
for both IPF (solid symbols) and OPF (open symbols) flames. Experimental microgravity results of current 
study are shown (♦). Current numerical results for IPFs and OPFs are shown as triangles (A and ▲ ). 
This set of data is compared to numerical data obtained for IPF and OPF flames of Sun et al. [1], shown 
as circles (o and • ). Note that their OPF data were modified by using our definition of burning velocity. 
Other experimental results for OPFs of Tseng at al. [3], Taylor [26] and Palm-Leis et al. [4] are also shown 
for comparison (0)- All available data show a definite difference between the Markstein values associated 
with OPFs and IPFs. 
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Figure 5; Propane- Air IPF Case (<£ — 0.6): Computed temperature and molar concentrations of choice 
species across the outwardly propagating flame. The reference unstretched case is shown as the dotted lines. 
The solid lines represent the flame structure for the stretched case (K = -704s -1 ). The x - 0 location 
corresponds to the point of maximum heat release. The origin is located to the left of the plot. Notice that 
stretch causes a moderate increase in final temperature and intermediate species such as H, OH and CO; 
this increase is linked to a large increase in burning velocity (because of the positive Ma). The increase in 
intermediate species is also associated with an increase in heat release. For the stretched case there is also a 
very noticeable increased concentration of H 2 ahead of the flame. The flame thickness, as defined from the 
temperature profile, is greatly diminished for the stretched flame. 
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Figure 6: Propane- Air OPF Case ($ = 0.6): Computed temperature and molar concentrations of choice 
species across the outwardly propagating flame. The reference unstretched case is shown as the dotted lines. 
The solid lines represent the flame structure for the stretched case (K = +113s~ 1 ). The x = 0 location 
corresponds to the point of maximum heat release. The origin is now located to the right of the plot. Notice 
that stretch causes a small decrease in final temperature and intermediate species; this decrease is linked to 
a decrease in burning velocity (because of the positive A/a). The decrease in intermediate species is also 
associated with an decrease in heat release. 
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Figure 7 : Methane- Air Results: Compilation plot showing numerical results obtained for both IPF (solid 
symbols) and OPF (open symbols) flames. Numerical results for complex chemistry IPFs and OPFs are 
shown as triangles (A and A ). Results obtained for simple step chemistry are shown as diamonds ( 0 and 
4 ). Note that IPF Markstein values are larger then OPF values. Experimental OPF data of Tseng et 
al. [3,25] are shown for comparison ( o ). 
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Figure 8: Lewis number dependence for simple chemistry ( Methane-air; $ = 0.66 ): Numerical 
results of Markstem numbers versus Lewis number of the fuel for methane/air flames using simple one- 
step chemi: stry. Both plots show results for IPFs (♦) and OPFs (0). Figure 8(a) shows the dependence 
of Markstem number on the fuel Lewis number for the simple mechanism. Notice there is still a large 
discrepancy between values obtained for IPFs and OPFs. Figure 8(b) shows results when heat release is 
owere to negligible values. In this case, both flame configurations are purely curved. Both Markstein values 
collapse onto each other for a wide range of Lewis numbers. 
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